Thomas-Fermi-Poisson theory of screening for laterally confined and unconfined 
two-dimensional electron systems in strong magnetic fields 
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We examine within the self-consistent Thomas-Fermi-Poisson approach the low-temperature 
screening properties of a two-dimensional electron gas (2DEG) subjected to strong perpendicu- 
lar magnetic fields. Numerical results for the unconfined 2DEG are compared with those for a 
simplified Hall bar geometry realized by two different confinement models. It is shown that in the 
strongly non-linear screening limit of zero temperature the total variation of the screened potential 
is related by simple analytical expressions to the amplitude of an applied harmonic modulation 
potential and to the strength of the magnetic field. 
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I. INTRODUCTION 

A two-dimensional electron gas (2DEG) in a 
strong perpendicular magnetic field has unusual low- 
temperature screening properties, 1,2 since the highly de- 
generate Landau-quantized energy levels lead to a strong 
variation of the thermodynamic density of states (TDOS) 
with varying strength of the magnetic field, i.e., with 
varying filling factor v of the Landau levels (LLs). If 
a LL is close to half filled, the TDOS is very high (in- 
versely proportional to the temperature T), and static 
potential fluctuations are nearly perfectly screened. We 
will consider only spin-degenerate 2DEGs, so that this 
happens if the value of v is close to an odd integer, while 
at even-integer v the Fermi energy lies in the gap be- 
tween two adjacent LLs and a spatial redistribution of 
electrons and, therefore, a screening of (weak) potential 
fluctuations is impossible. In an inhomogeneous 2DEG 
with sufficiently strong long-range density fluctuations, 
screening effects lead to quasi metallic (so called "com- 
pressible") regions with high TDOS, in which screening 
is nearly perfect and a LL is "pinned" to the Fermi 
energy, and to insulator-like "incompressible" regions, 
which separate adjacent compressible regions. In the in- 
compressible regions the Fermi energy falls into the gap 
between two LLs and the electron density n e i(r) is con- 
stant (even- integer filling factor), while in the compress- 
ible regions n c \(r) adjusts itself so that the self-consistent 
electrostatic potential energy V(r) of an electron differs 
from the Fermi energy (more precisely the electrochemi- 
cal potential /x*) by a Landau energy hoj c (n + l/2), where 
lu c = eB/m is the cyclotron frequency in the magnetic 
field B. As a consequence, V(r) becomes nearly constant 
within a compressible region and differs by integer multi- 
ples of hio c between different compressible regions. Lan- 
dau level pinning and the interplay of compressible and 
incompressible regions lead to strongly nonlinear screen- 
ing effects. This screening scenario has been established 
some time agoi* and was applied, e.g., to calculate, at 
zero temperature, the electronic DOS* and transport** 
through 2DEGs in smooth periodic and random poten- 
tials. The explanation**^ of several experimental results, 



e.g., on quantum Hall devices under high currents close 
to the breakdown of the quantized Hall effect j*i* rely also 
on these ideas. A systematic investigation of these inter- 
esting nonlinear screening effects is, however, apparently 
not available in the literature. 

Models for half-space and Hall-bar geometries with 
planar charge distributions have been proposed that al- 
low closed solutions of Poisson's equation (i.e., the cal- 
culation of the potential for given electron density), and 
estimates of position and widths of the incompressible 
strips have been given^*** By adding the non-linear 
Thomas-Fermi approximation for the calculation of the 
electron density from the potential, that work was ex- 
tended to a self-consistent approach, which allows us to 
calculate both electron density and electrostatic potential 
for arbitrary temperature^*** This approach, which we 
will employ in the following, shows that the existence and 
the width of incompressible strips depends sensitively on 
temperature, and allows to calculate their position and 
width for given background charges without additional 
assumptions. 

The purpose of the present work is a systematic inves- 
tigation of the nonlinear low-temperature screening of 
harmonic electrostatic potential modulations in laterally 
confined and unconfined 2DEGs subjected to a quantiz- 
ing perpendicular magnetic field. We will demonstrate 
that in general edge effects do not qualitatively change 
the screening properties of the 2DEG, even if the sample 
width is not much larger than the period of the imposed 
potential modulation. There are, however, peculiar dif- 
ferences between confined and unconfined 2DEGs in sit- 
uations, in which the latter have no states near the Fermi 
energy. To understand this in detail, we first discuss the 
screening of a potential modulation imposed on a homo- 
geneous 2DEG (Sec. II) and then consider, for two differ- 
ent boundary models, edge effects on screening in Hall 
bar geometries (Sec. III). 

We will assume the 2DEG to be located in the plane 
z = with a (surface) number density n e \(x) and con- 
sider only situations with translation invariance in the y 
direction. The (Hartree) contribution Vh(x) to the po- 
tential energy of an electron caused by the total charge 
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density of the 2DEG can be written &ai& 
2e 2 



V H (x) 



dx'K(x, x') n \(x'), 



(1) 



where — e is the electron charge, R an average back- 
ground dielectric constant^ and the kernel K(x,x') de- 
scribes the solution of Poisson's equation with appropri- 
ate boundary conditions at xi and x r . The electron den- 
sity in turn is calculated in the Thomas-Fermi approxi- 
mation (TFA)A2 



{x) = J dE D(E)f([E + V(x) - n*]/k B T) 



(2) 



with D(E) the relevant (single-particle) density of states 
(DOS), /(e) = [1 + e e }~ 1 the Fermi function, /z* the elec- 
trochemical potential, and with V(x) = V ex t(%) + Vh{%) 
the total potential energy of an electron, which differs 
from Vh(x) by the contribution due to external charges, 
e.g., a homogeneous positively charged background and 
a charge distribution creating a periodic modulation po- 
tential. The local (but nonlinear) TFA is much simpler 
than the corresponding quantum mechanical calculation 
and expected to yield essentially the same results if V(x) 
varies slowly in space, i.e., on a length scale much larger 
than typical quantum lengths such as the extent of wave 
functions or the Fermi wavelength. 



II. 



HOMOGENEOUS 2DEG 



We start with a homogeneous 2DEG described by the 
DOS D (E) = DqO(E), with D a = m/{-nh 2 ), for B = 0, 
and by the Landau DOS 



D B {E) 



l^E-EJ, E n 

m n=0 



hio c (n + 1/2) (3) 



for finite B. The effective mass of an electron is de- 
noted by m, the magnetic length by l m = yjh/ (mw c ). In 
both cases we assume spin degeneracy and neglect col- 
lision broadening effects. The constant electron density 
n c \ is given by Eq. J5J with n c \ = n c \(x), V(x) = and 
/j* = /i the chemical potential. For these simple models 
the energy integral in Eq. J5J is readily carried out. We 
tacitly assume that the electron charges are neutralized 
by a homogeneous background of positive charges. 



A. Kernel for periodic modulation 

We now add a periodic external modulation described 
by a potential energy V ex t(%) = V ex t(x + a). The 2DEG 
will respond with a density modulation and a Hartree 
potential of the same period a. To exploit the periodic- 
ity, we expand density and potentials into Fourier series 
according to 

^ C a l 2 rW 

V{x)^y^V q e iqx 1 V q = —e- iqx V{x), (4) 

J -a/2 a 



with q — 2irn/a and integer n. To maintain charge neu- 
trality we require — n e i in any case. With the bound- 
ary conditions Vh{x, z) — > for \z\ — > oo, Poisson's equa- 
tion yields (see, e.g., Ref. 0) 



V q (z) = (2ne 2 /R\q\)e-^n q cl 



(5) 



as response to the density fluctuation n q r Summing over 
harmonics (for q ^ 0)^ we obtain Vh(x, z — 0) from 
Eq. JQ) with — xi = x r = a/2 and the kernel 



K(x,x') = -In 2 sin -fx - x 1 ) . (6) 
a 



B. Breakdown of linear screening 

1. Zero magnetic field 

In the limit B = 0, T -> and with E F = ^*(B = 
Q,T= 0), Eq. © reduces to 



n d {x) = D a (E F - V{x)) 9(E F - V(x)), 



which is a linear relation between V(x) and n e \{x) for 
V(x) < Ep. With Eq. (J5J we find for a harmonic poten- 
tial modulation V ext {%) = V q Kt cos qx a harmonic density 
modulation Sn e \(x) = n^cosqx and the self-consistent 
("screened") potential V(x) = I^^cosga; with 



V q 



e(q) = l + Q /\q\. 



(8) 



The dielectric function e(q) can be expressed in terms 
of the effective Bohr radius ag = Kh 2 /(me 2 ) (for GaAs 
a* = 9.8 nm), since Q = 2ne 2 D /k = 2/a^ii& With 
q = 2Tr/a, the screening strength is thus determined by 
the dimcnsionless parameter 



= na^/a. 



(9) 



We will assume a -C 1, i.e., e(q) = l + l/a> 1, so that 
the TFA is valid for B > 1 T, i.e. l m < 30 nm. Since 
in the linear screening regime the minimum value of the 
electron density is n e i(0) = D (Ep — V q ), linear screen- 
ing breaks down if the modulation strength becomes so 
large that V q > E Fl i.e., for V q xt > e(q)E F . For larger 
modulation amplitude the redistribution of electrons is 
hindered: while the electron density at the minimum of 
V cxt (x) still increases, the electron density at the max- 
imum of Vext(x) cannot decrease further. Instead the 
density minimum becomes broader. This means that the 
electrons are depleted from strips along the maxima of 
V ex t(x) and the 2DEG breaks off into a system of par- 
allel quasi one-dimensional ribbons. Thus, the imposed 
harmonic modulation potential V cx t(x) now leads to an 
anharmonic density distribution and, therefore, an an- 
harmonic screened potential. 17 Mathematically, Eq. (JTJ 
with 10 and Eq. J7J now represents a nonlinear integral 
equation that must be solved numerically. In Fig. ^ w e 
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FIG. 1: Variance of the screened potential versus ampli- 
tude of the harmonic potential modulation imposed on a 
spin-degenerate homogeneous 2DEG with a half filled Lan- 
dau level, for several odd-integer values of the filling factor. 
Default temperature k B T/E F = 0.001, e(q) = 41. 



plot the total variance Var[V] = V(Q) — V(a/2) as a func- 
tion of the amplitude V^ xt = Vq of the imposed modula- 
tion potential for several values of the magnetic field. The 
result for B — and T = is shown as a thick solid line. 
In the linear screening regime, Var[V] = 2V q = 2Vo/e(q). 
As linear screening breaks down, a kink appears in the 
line and the variance increases much faster than in the 
linear regime. With increasing temperature this kink is 
rounded off, while the Var[V]-vs-Vb curve as a whole is 
not much affected (shown for k B T/Ep = 0.04 by the 
open circles in Fig. [I}. Here and in the following we mea- 
sure energies in units of the Fermi energy Ep = n B \/ Dq 
(for GaAs with n cl w 3 • 10 11 cm~ 2 , E F rj lOmeV), and 
we keep the mean electron density n c \, and thus Ep, 
constant. We will focus in the following on the regime 
Vq < e(q)Ep, where screening is linear in the limit B = 0, 
T = 0. 



2. Half filled Landau levels 

With the Landau DOS [see Eq. ©] and the definition 
of a position-dependent chemical potential, fJ,(x) = /i* — 
V^x), Eq. © yields 



n cl (x) = huj c D Q ^ f( E n ~ K x )) ■ 



v n » 2f(E n — fi) is the filling factor of the n-th Landau 
level, the TDOS is± 



D T (fi;B) 



dn. 



djjL 



k B T 2 



(i-fK (") 



which is peaked around fi = E n with a maximum value 
Huj c Dq/ (Ak B T) and a width of order fcgT [a crude ap- 
proximation is D T (n;B) w (hu) c D /4k B T)6(2k B T - 
\E n — fi\)]. Linearizing Eq. (|10|) with respect to the 
screened potential V(x), we obtain the Eq. © with Qq 
replaced by Q B = QqD t {^; B)/D > Qo, 



B) = l + —— y (i - yj tt 



(12) 



For exactly half filling, /i = E n , his is a rather good ap- 
proximation, as can be seen from Fig. ^ which shows 
numerical results for v = uq = 1 (n = 0) at two different 
temperatures (two lowest curves at Vo/Ep > 15). We 
see that the screening at finite magnetic field depends 
much stronger on temperature than at B = 0. The lin- 
ear approximation breaks down, if the amplitude of the 
screened potential becomes of the order 2k B T. 

This yields, in the limit of low temperatures and for 
v n = 1, the estimate for the linear screening regime, 



Hip hip v 



(13) 



(10) 



with e(q) = e(q; B — 0) = l + l/a. For a larger modula- 
tion the redistribution of electrons within the considered 
LL is not efficient enough to screen the imposed modula- 
tion potential, and similar to the B = case, the variance 
of the screened potential increases much stronger than in 
the linear regime. In Fig. ^ we show low-temperature 
(k B T/Ep — 0.001) results for odd-integer v values, cal- 
culated numerically from Eqs. 0J, (JHJ), and (|10|l . For this 
temperature, the linear increase of the screened potential 
with the applied modulation amplitude Vo is not resolved 
on the scale of Fig. ^ However the rapid increase of the 
variance of the screened potential at Vo/Ep ~ e(q)/P is 
clearly seen for the indicated v values. 

For v = 1 the situation is very similar to the B = 
case, apart from the fact that screening in the linear 
regime is much stronger ( "perfect screening" , "pinning of 
lowest LL to Fermi level") due to the higher DOS. For 
v > 1 new phenomena occur, which we will now discuss. 



C. Emergence of incompressible strips 



We may also write the argument of the Fermi function 
as E„(x) — fi* and interpret E n (x) = E n + V(x) as 
position-dependent Landau energies, which is correct if 
the Thomas-Fermi approximation holds. It will be use- 
ful to define, in addition to the average LL filling factor 
v = 27r£^ l n e i, a local filling factor v(x) = 2~Kl^ n n e \{x). 

For k B T <C hu> c and a homogeneous 2DEG with partly 
filled n-th Landau level, \x ~ E n , v = 2n 



v n , where 



1. Odd-integer filling factor v 

We start with filling factor v = 3 and investigate the 
changes of the electron density (Fig. |2J) and of the total 
potential (Fig. 0) with increasing amplitude Vq of the 
imposed modulation V ex t(x) = Vbcosgx, giving explicit 
results for the six typical Vq values indicated by open 
circles in Fig. [21 
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FIG. 2: Variance of the total potential versus Vb, for average 
filling factor v — 3. The insets show the local filling factor 
v{x) in one modulation period (0 < x/a < 1) for the six 
Vb values indicated by circles. Parameters: ksT/EF = 0.01, 
e(q) = 41, q = 2ix/a. 
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FIG. 3: Total potential (thick solid lines) and the three lowest 
of the corresponding Landau levels (thin solid lines) together 
with the electrochemical potential (thick dashes lines) for the 
six Vb values indicated in Fig. [5] Parameters as in that figure. 



For Vb = 11 Ef, close to the breakdown of linear 
screening (case 1), the density is strongly modulated (see 
thin line in lower inset of Fig. [2Jl, but the modulation 
appears still cosine-like. The potential is so effectively 
screened that the second- lowest LL (n = 1) is pinned 
(within a few fcgT) to the Fermi energy [see Fig. Efl)]- 
For case 2, Vb = 13.7 Ep, the total potential has de- 
veloped locally confined maxima and minima, while it 
remains rather flat in between [Fig. Et2)]. Near these ex- 
trema \Ei(x) — /i* | becomes so large that the LL n = 1 



is completely occupied (near x = a/2) or empty (near 
x = 0), and incompressible strips with local filling factors 
v{x) = 4 and v[x) — 2 develop near the potential minima 
and maxima, respectively [see thick solid line in the lower 
inset of Fig. [2 and Fig.[3f2)]. Increasing the modulation 
to Vb = 14.5 Ep (case 3) leads to more pronounced local 
extrema and broader incompressible strips, but does not 
change the situation qualitatively. The overall change of 
the density distribution is rather small (see lower inset of 
Fig. indicating poor screening. Indeed the slope of the 
Var[V]-vs-Vb curve in this regime is AVar/AVb w 1, i.e., 
only slightly smaller than in the absence of any screen- 
ing, which would yield AVar/AVb = 2. We note that 
in the incompressible strips the local filling factor (i.e., 
the density) is constant, while in the pinning regions, i.e., 
the compressible strips, the screened potential still has a 
finite slope, proportional to ksTm^ 

As Vo increases further to case 4 (Vb = 17 Ep), the 
modulation becomes so strong that the maximum of the 
lowest LL, Eo(0), and the minimum of the lowest unoc- 
cupied LL, E2(a/2), reach the Fermi level fi* (to within 
ksT). Then thermal population of the higher LL (n = 2 
near x = a/2) and depletion of the lower LL (n = 
near x = 0) starts and compressible strips emerge in 
the center of each incompressible strip [see dashed line 
in the upper inset of Fig. [21 and Fig. Ef4)]. Further in- 
crease of Vb up to Vb = 38.5 Ep (case 5) widens the 
compressible strips and leads to a strong increase of the 
density modulation due to a redistribution of electrons 
from the n = to the n = 2 LL. This results in a 
strong screening, similar to that in the linear screening 
regime at weak modulation, and, apart from a weak in- 
crease with a slope proportional to ksT, the variance 
Var[V] = V(0) — V(a/2) remains constant at the value 
Var[V] = fi* — E o (0) - [fi* - E 2 (a/2)} = 2hu c . This 
plateau behavior of the Var[V]-vs-Vb curve is obviously 
an immediate consequence of the pinning of LLs to the 
Fermi level, i.e., of the nearly perfect screening. 

In case 5 we reach a situation in which the lowest 
(n = 0) LL is nearly empty at the potential maximum 
and the higher (n = 2) LL nearly full at the potential 
minimum. In case 6 (Vb = 41.2 Ep) the total potential 
again develops local extrema, similar to the situation de- 
picted in Fig. 0^2). But now the incompressible strip 
created at the potential maximum is due to the depop- 
ulation of the lowest LL, i.e., due to vanishing electron 
density. With further increasing Vb the depletion regions 
become wider and the density near the potential minima 
increases, but screening remains much poorer than in the 
plateau region. 

We see from Fig. [21 that the global appearance of the 
density modulation, apart from a fine-structure related to 
the incompressible strips, is more or less cosine-like. We 
will use this finding for a rough estimate of the plateau 
width of the Var[V]-vs-Vb curves. First we conclude from 
the cosine-like form of the induced density variation that 
in the high-screening plateau region, along with Eqs. @ 
and (jHJ , the first relation of Eq. (|13|) holds qualitatively 
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and relates the changes SV of the total potential to the 
changes SV of the externally applied potential by SV ~ 
e(q;B)SVo. For Q/\q\ S> 1, this yields for the change of 
the variance Var[V] across the plateau of width AVq 



AVar[V] 



8k B T 
hw c e(q) 



AVq 



(14) 



i.e., an estimate for the slope of the Var[V]-vs-Vo curve 
in the plateau region. Since the modulation induces den- 
sity changes Sn e \ mainly within the compressible regions 
of high TDOS, we estimate Sn c \ ~ —Dt{h; B)8V ~ 
-(hu) c /4Je B T)D SV (which holds for u n = 1 and \SV\ < 
2k B T). In terms of Sv = 2irl'^ n 8n e \ this yields SV ~ 
2k B T5v : and together with Eq. (|14|) the relation 



£ F 2P L J 



(15) 



between the plateau width AVo and the change of the 
filling factor variance Var[j/] = v(a/2) — v{0) (defined at 
fixed Vq) across the plateau. This criterion applies also to 
the small- Vq linear-screening regime, in which v(x) varies 
within the same LL, with Var[^] increasing from to 2, 
i.e. AVar[z/] = 2 [sec Eq. l(T3|)], The resulting width of 
the linear-screening regime, AVq/Ef ~ e(q)/D, describes 
the numerical results of Fig.nf or odd-integer v = 2n + 1 
quite well. From the discussion of Figs.[21and|31we expect 
that, for n > the linear regime of the Var[V](Vb) curve 
is terminated by a step of height 2hw c = AEp/i>, which 
is followed by a plateau. While Vq sweeps through the 
plateau, in addition to the LL with index n the two LLs 
with indices n — 1 and n + 1 are locally pinned to the 
Fermi level and lead to a total change AVar[^] = 4 [for 
v = 3 from Var[^] = v(a/2) — v(0) — 2 on the left side 
to Var[t/] = 6 on the right side of the plateau, as seen 
from the upper inset of Fig. |2J. This yields the plateau 
width AVq/Ef w 2e{q)/v. If n — 1 = 0, the plateau will 
be followed by the breakdown regime. If n — 1 > 0, the 
plateau will be followed by a further step of the same 
height to a plateau of the same width. 

To summarize: at very low temperatures and odd- 
integer filling factors v = 2n + 1, the variance Var[V] of 
the screened potential as function of the imposed modu- 
lation amplitude Vq shows a linear screening regime for 
Vo < ^{q)Epjv which is followed by n successive steps 
of height 2hw c = 4E F /D and width AV w 2e(q)E F /D. 
The plateau of the n-th step ends at the breakdown 
of the 2DEG into a pattern of isolated ID systems, 
which leads to poor screening and is indicated in the 
Var[V / ](Vb) curve by a slope of order unity. At finite 
temperature, the plateaus assume finite positive slopes, 
which are estimated from Eq. I|14l) as AVar[V]/AVo ~ 
Avk B T '/[e(q)EF\- These results, which describe the con- 
tent of the numerically calculated Figs. lll3l verv well, de- 
pend, of course, on the high symmetry of the situations 
considered so far. 



2. Even-integer filling factor v 

Another situation of high symmetry is that of an even- 
integer filling factor v = 2n + 2, where n is the in- 
dex of the highest occupied LL and the Fermi energy 
Ep = huj c (n + 1) lies in the middle between two adja- 
cent LLs. According to Eqs. (|13l) and (|12fl . the linear 
screening regime shrinks to zero, since v n — 2. Thus, 
at very low temperature (fc^T hw c /2) a weak modu- 
lation V cx t(x) — Vocosqx will not be screened, i.e., the 
local filling factor will be independent of the modula- 
tion, v{x) = 9, and the total potential will equal the 
external one, with variance Var[V / ](Vb) — 2Vq. This sit- 
uation changes when the modulation potential becomes 
so large that \Vq — fku c /2\ ~ /csT, i.e., the maximum en- 
ergy -EVi(O) of the highest occupied, and the minimum 
energy E n+ \{a/2) of the lowest unoccupied LL approach 
the Fermi energy. Then, with increasing Vq the LL n is 
depleted near x = while the LL n + 1 is populated near 
x = a/2, forming compressible strips with local filling 
factors v{x) < v and v(x) > D, respectively. 

This is demonstrated in the inset of Fig. 0] which shows 
the local filling factor v(x) for ksT/Ep = 0.001 and the 
average filling factor D = 2, i.e., fku c = Ep, and for the 
modulation strengths Vq indicated by circles in the main 
figure. For Vq/Ef = 0.45 in the non-screening region 
the deviation [^(x) — v\ is practically zero (numerically 
< 10~ 6 ), while for V /E F = 0.6 it is finite, although 
small (in the inset enhanced by a factor of 20), with 
narrow compressible strips. Between Vq/Ef = 0.5 and 
Vq/Ef ~ 41 the width of the compressible strips and the 
deviation [v{x) — D] increase continuously, while the vari- 
ance Var[y] sa ftw c — Ep remains constant. Since screen- 
ing is due to the redistribution of electrons at the Fermi 
energy, i.e., to electrons in the compressible strips where 
the TDOS is large, we may again use Eq. (|15fl to estimate 
the plateau width. At the beginning of the first plateau 
the filling factor is constant, v{x) = D, i.e., Var[^] = 0. 
At the end of the plateau, the LL n is depleted at the po- 
tential maximum, i/(0) = D — 2, and at the potential mini- 
mum the LL n + 1 is full, v(a/2) — D + 2, i.e., Var[i/] = 4. 
Thus, we have to use Eq. (|15H with AVar[j/] = 4 and 
obtain for the plateau width AVo ~ 2e(q)Ep/i'- 

This estimate is obviously in good agreement with the 
numerical calculations presented in Fig.0] For filling fac- 
tor v — 2, i.e., n = 0, the first plateau ends at the transi- 
tion to the poor-screening breakdown regime, since then 
the lowest LL n = is completely depleted at the poten- 
tial maxima (thick dash-dotted line in inset of Fig. 0J. 
For v — 2n + 2 with n > a behavior similar to that dis- 
cussed in Fig. |3 occurs. As Vq increases slightly beyond 
the plateau regime, a narrow local maximum of V{x) de- 
velops near x = 0, accompanied with an incompressible 
strip of filling v(x) — 2n due to the local depletion of the 
LL n. Simultaneously, a narrow local minimum of V(x) 
develops near x = a/2, accompanied with an incompress- 
ible strip of filling v(x) = 2n + 4 due to the local occupa- 
tion of the LL n + 1. Then, in a narrow Vq interval these 



6 



4 F 



-J— 


0.45 

0.60 *) 

10.0 

35.0 

41.2 


- / \ - 

/ \ 






10 



20 

V„/E F 



30 



40 



3.0 



2.0 



■ft 
> 

o 

> 



1.0 - 



0.0 



-- - v= 4/3, £2=1.5 

v= 7/3, £2=0.857 

v= 8/3, £2=0.75 

v=13/4, £2=0.615 

v=15/4, £2=0.533 

v=15/2, £2=0.267 



---}- 
—fT~ 



li 

l! 



10 



20 

V„/E c 



30 40 



FIG. 4: As in Fig. but for even-integer values of the (av- 
erage) filling factor (k B T/E F = 0.001, e(q) = 41). The thin 
solid line indicates the result for B — 0, T = 0, the thin 
dashed line has slope 2. The inset shows the local filling fac- 
tor v(x) in one modulation period for average filling v = 2 
and the five values of Vo indicated by circles in the main fig- 
ure. *) For Vo/Ef = 0.60 the deviation [u(x) — 2] is enhanced 
by a factor 20. 



new extrema become more pronounced and the accom- 
panied incompressible strips widen a little. However, the 
accompanied density change is small, resulting in a poor 
screening and a rapid increase of the Var[V](Vb) curve. 
This interval ends when the new maximum E n -i(0) of 
the Landau level n — \ and the new minimum E n+ 2{a/2) 
of the LL n + 2 come close to the Fermi level (within a few 
ksT). Then, with further increasing Vo, new compress- 
ible strips open at the locations of the potential extrema, 
and a plateau region of the Var[V](Vb) curve with "per- 
fect" screening sets in. We thus again find a step behavior 
like in Fig.[T]with step height AVar[V~] = 2hw c = 4E F /D. 
During the Vq sweep through the corresponding plateau 
the LL n—l will be depleted near x$ while the LL n + 2 
is occupied near x = a/2. Thus, we can estimate the 
plateau width from Eq. (jPo|> with AVar[i/] = 4. The last 
plateau is the one corresponding to the local depletion of 
the n = LL. 



In summary, for v = 2n + 2 and very low temperature, 
the Var[V](Vo) curve shows a linear increase with slope 
2 for < Vb < Ep/P, followed by a plateau of height 
hoj c = 2Ep /v and width AVb ~ 2e(q)Ep/D. This plateau 
is followed by n steps of height 2fttu c and approximately 
the same width AVb- The plateau of the last step is 
followed by the breakdown regime. 



FIG. 5: As in Fig. but for some non-integer values 
of the (average) filling factor and for higher temperature, 
k B T/E F = 0.01 ( e(q) = 41). 



3. Non-integer filling factor v 

In Fig. [SI we show Var[V](Vo) curves for a few non- 
integer values of the average filling factor v = 2n + v ni 
with < v n < 2. Although these results may, at a first 
glance, look confusing, we will now demonstrate, that 
they can easily be understood, and even predicted, from 
a few simple principles. 

To estimate the width of the linear screening regime 
at small Vq values, we follow the reasoning of Sec. Ill Bl 
V^t ~ e i a \B)V q 1 but we note that the linear ap- 
proximation to the Taylor expansion of n c \(x) with re- 
spect to V{x) [see Eq. (|TU)l ] does no longer hold for 
|V(x)| ~ 2ksT, since for v n ^ 1 the second order term 
[oc d 2 n c i/dfi 2 = (dn c i/d[i)(l — v n )/kBT] yields already 
noticeable contributions for smaller V(x). To take this 
into account, we use the linear approximation only for 
| V(x)| < 2ksT/ (1 + |1 — f n \) and obtain as condition for 
the linear screening regime 

Vo_ < e(g) v n {2 - v n ) 

E F ~ v l + |l-i/„| ' [ ' 

For v n — 1 this reduces to the estimate (|13|l ■ But in addi- 
tion, Eq. I |ltj| l states that for even-integer filling, v n — > 
or v n 2, the linear screening regime shrinks to zero, 
and it provides a good description of the widths of the 
linear screening regimes in the examples shown in Fig. [HI 
We will now use Eq. H15|) to obtain estimates of the 
plateau widths and heights of the Var[V]-vs-Vb curves, 
which contain the estimate IjlrJI) for the linear screen- 
ing regime as a special case. Our estimates are based 
on the observation that in all cases we have studied 
the density modulation is nearly symmetric about the 
average density, so that the average of extreme values 
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of the local filling factor is close to the average filling, 
u(0) + v(a/2) w 2D. Nearly perfect screening occurs, 
if at both the potential maxima and the minima a LL 
is pinned to the Fermi level, so that electrons can eas- 
ily be redistributed between these LLs and both u(Q) 
and v(a/2) arc different from even integers. If, with 
increasing Vo, ^(0) approaches an even integer value, 
v(Q) = 2k, the LL k is depleted at x = and a local 
maximum of V(x) starts to develop there. Screening re- 
mains poor, and Var[V](Vo) increases rapidly, until the 
LL k — 1 reaches the Fermi level at x = 0. Then a step 
of height hio c is completed and the next plateau with 
perfect screening starts. A similar step begins as v(a/2) 
reaches the value 2k'. Then the LL k 1 — 1 is completely 
filled at x — a/2 and a local potential minimum starts 
to develop there. Perfect screening begins again, if the 
LL k' reaches Fermi level at x = a/2 and the increase of 
Var[V](Vb) by fku c is completed. 

We combine now these consideration with the estimate 
(|f 5fl . For convenience, we introduce the dimcnsionlcss 
variables 



V 



Vo 



e(q)E f 



Ej 



(17) 



and focus on the regime < vo < 1, in which for T = 
and B = screening is linear and leads to Var[u](t>o) = 
2vo. To keep the discussion simple, we consider the two 
possible cases of non-integer D — 2n + v n separately. 

a. For < v n < 1, the end of the linear screening 
region, where v <C f2, is reached when v(0) — 2n. Then 
v(a/2) k, v + v n and Var[V] = 2v n , and across the linear 
screening regime we find AVar[^] = 2v n . According to 
Eq. H15JI . linear screening ends at Vq/Ef ~ e(q)v n /D, in 
agreement with Eq. (|16[1 . Neglecting potential variations 
oc kWT <C Q, we note 



Var[u] « , if < v < 1 - ntt , 



(18) 



since v n = D — 2n and il — 2/ v. If n = 0, larger Vq lead 
to the poor-screening quasi ID ribbon regime. 

For n > 0, the next plateau terminates when v(a/2) = 
2n+2. Thcn^(O) « 2n-2(l-v n ), i.e., Var[z/] = 2(2-v n ). 
Across this plateau we have AVar[j/] — 4(1 — u n ), and 
with Eq. (CSJ) AV a w 2e(q)E F (l - v n )/D. This yields 
Av Q w (2n +l)n~2 and 

Var[w] w fi , if 1 - nfi < u < (ra + l)fi - 1 . (19) 

This plateau is followed by another one along which v{§) 
decreases to 2(n — 1), while v{a/2) increases to rj 2n + 2 + 
2;y„ and thus Var[i/] to A + 2v n . Thus, across that plateau 
we find AVar[^] = Av n and A?;o = 2 — 2ni7, which leads 
to 

Var[w]w2rj, if (n+l)fi-l< u < l-(n-l)fi . (20) 

If n = 1, this plateau is followed by the poor screening 
quasi ID ribbon regime. If n > 1, we are in the same sit- 
uation as at the end of the Iow-Vq linear screening regime, 



and a double step of total width f2, consisting of one step 
of height fl and plateau width Avq = (2n + l)fi — 2 and 
another one of height and width Avq = 2 — 2nfl, will 
follow. Thus, we obtain for < k < n 

Var[u] w (2k - l)fi , if (21) 
1 - (n + 1 - fc)n < w < (n + fc)f2 - 1 , 

Var[w] w 2fcO, if (22) 
(n + k)Q -l<Vo<l — (n— k)Q . 

Thus, the linear screening regime is followed by n double 
steps, which sum up to a total width Avq — 1, and on 
the last plateau (before breakdown) we have Var[w] = 
2(1 — v n /v). For < v n < 0.5, the first plateau of the 
double step is wider than the second one, as for the dash- 
dotted line in Fig. [S] (y = 2.33), while for 0.5 < v n < 1 
the second plateau of the double step is the wider one, 
as for the short-dashed line [y — 2.67). 

We should mention three limits. For v n — > the low- 
Vq linear screening regime shrinks to zero and the first 
plateau of the double step exhausts its full width, so that 
the second step merges with the first one of the following 
double step. Thus, we observe at small Vq a step of height 
Q, followed by steps of the double height 2S1, and all 
plateaus have the same widths, as we found previously. 
For v n — 0.5 we get an even number of steps which all 
have the same heights and widths. For v n — > 1 the width 
of the first plateau of each double step shrinks to zero, 
so that the low-Vo linear screening regime is followed by 
steps of height 2il and width J7, as we have seen before. 

b. For 1 < v n < 2, we have at the end of the linear 
screening regime v(a/2) = 2n + 2 and i/(0) ks 2n + 2(f n — 
1) > 2n, with Var[t/] = 2(2 - v n ). From Eq. JT^ we 
obtain 



Var[w] »0, if < v Q < (n + l)fi - 1 . 



(23) 



and see that now the linear screening regime is always 
followed by a step of height £1 to a plateau of perfect 
screening. For n = (i.e., 1 < O < 2) this plateau covers 
the interval f2 — 1 < vq < 1 (see e.g., thin dashed line 
of Fig. [S] for v — 1.33). To estimate for n > height 
and width of the following steps and plateaus, respec- 
tively, we proceed as before, exploiting that at the end 
of each plateau either ^(0) or v(a/2) reaches an even in- 
teger value, and that v(0) + u(a/2) w 2v. The result is, 
for < k < n, a double step of total width fl, 

Var[t>] « (2k - if (24) 

(n + k)Cl - l<u <l-(n. + l- fc)0 , 

Var[u] « 2kQ, if (25) 
1 - (n + 1 - k)n < v Q < (n + k + l)fi - 1 , 

which is followed by a final single step, 

Var[u] « (2n + , if (2n + 1)0 - 1< v < 1 , (26) 

of height fl and plateau width 2 — (2n + 1)0. The linear 
screening regime [width (n+ l)f2 — 1], the n double steps 
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and this final plateau cover together the interval < vq < 
1, as in the case < v n < 1. The variance of the screened 
potential in the last plateau is Var[i;] = 2[1 — (v n — f ) jv\. 

For 1 < v n < 1 .5 the first plateau of each double step 
is wider than the other plateaus, as seen in Fig.[5]for the 
long-dashed line (v = 3.25), which exhibits one double 
step following the initial single step. For 1.5 < v n < 2 
these first plateaus are the narrower ones, as seen for the 
thick solid line (p = 3.75). For v n = 1.5, the double 
steps consist of two individual steps of equal heights and 
plateau widths, as is illustrated by the thin solid line 
(P = 7.5)._ 

In the limit of odd-integer v, v n — > 1, the width of the 
first single step together with the second plateau width 
of each double step shrinks to zero, so that only steps 
with step height 2S1 and plateau width SI occur. (For 
v = 1, i.e., n = 0, no double step exists and the single 
step merges with the breakdown regime.) In the even- 
v limit, v n — > 2, the width of the first plateau of each 
double step shrinks to zero. Thus, the first step of width 
SI and height SI is followed by n steps of the same plateau 
width but double step height. 
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FIG. 6: The solid lines indicate Do = kQ — 1 and w = 1 — kfl 
for w = Vb/I-EpeCg)], Q = hu c /E F , and fc = 1, 2, . . . , 15. 
Even (odd) integer values of the average filling factor are in- 
dicated by dash-dotted (dashed) vertical lines. Within each 
area defined by the solid lines the value of Var [V}/Ef equals 
an integer multiple of SI. This value increases by fl, if a solid 
line is crossed in upward direction. The region vo > 1 corre- 
sponds to the poor screening regime of parallel, disconnected 
quasi ID electron systems. 

c. Summarizing the estimates of this Sec. Ill CI we 
note that the Eqs. (|18|l - (j26(l define a set of straight 
lines in the vq-Q plane, which separate areas in each of 
which the variance Var[u](i>o) equals an integer multiple 
of J7. This is schematically shown in Fig. El Position 
and height of the steps of the Var[w] (vq) curve for a given 
value of v can immediately be read off from this figure 



along the vertical line at SI = 2/ v. 

D. Sweeping the magnetic field 

We now consider the screening of an external cosine 
potential V ex t (%) — Vq cos qx of fixed amplitude Vq as 
function of the magnetic field B, keeping the average elec- 
tron density at the fixed value of the positive background 
charge density. Then, with increasing B the average fill- 
ing factor v = 2Ep/hu> c decreases. For the unmodulated 
2DEG (Vo = 0), this leads to the well known saw-tooth 
behavior of the chemical potential, which at low tem- 
peratures is pinned to the LLs, i.e. follows half-integer 
multiples of the cyclotron energy, 

li* = haj c (n+ 1/2) if l/(n+ 1) < SI < 1/n. (27) 

In the modulated 2DEG (Vo > 0), pinning of LLs to the 
electrochemical potential causes the total variance of the 
screened potential to be an integer multiple of the cy- 
clotron energy. Thus, for the variance Var[V] as function 
of B we expect a similar saw-tooth behavior as for the 
fi*-vs-B curve. Numerical results for several values of Vo 
are shown in Fig.[7| The uppermost curve for the largest 



2.0 




FIG. 7: Variance of the screened potential versus Sl c = 
hujc/Ep for Vo/Ef = 1, 5, 10, 15, 25, 35 (from bottom to 
top). The straight dashed lines indicate integer multiples of 
the cyclotron energy. (kBT/EF = 0.01, e(q) = 41). 

modulation amplitude looks indeed similar to a fi*-ys- 
B curve. However, whereas the latter with decreasing 
B always jumps to the next higher LL, the Var[V]-vs- 
B curves can also jump back to the next lower LL, as 
is more clearly seen for the curves with smaller mod- 
ulation amplitudes. This seemingly irregular behavior 
of the Var[V](Sl) curves in Fig. [7] can easily be under- 
stood from Fig. where we now have to follow hori- 
zontal lines. For fixed Vq = Vo/[Epe(q)] and decreasing 
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Q, the variance Vai[V)/ Ep increases by f2, if the hor- 
izontal line vq = Vq intersects one of the straight lines 
vq = kfl — 1, and decreases by il, if it intersects one of 
the straight lines vq = 1 — kCl. For large Vb (vq > 0.6), 
the variance jumps with decreasing B monotonically to 
higher multiples of the cyclotron energy, until B becomes 
so small, that il = 1 — «§. Then, for smaller B also 
jumps back to lower multiples occur. For small modu- 
lation amplitude (vq < 0.1) one has perfect screening if 
the average filling factor is not too close to an even inte- 
ger. Near such values linear screening breaks down and 
V&t[V]/Ef approaches fi, provided 2Vo > hw c (other- 
wise Var[V] = 2Vb). For sufficiently small B, of course 
the variance will equal higher multiples of the cyclotron 
energy, so that the correct linear screening limit is ob- 
tained in the limit of zero magnetic field. Thus, if one 
adds the smoothening effect of finite temperature, one 
can understand all properties of the apparently irregular 
Var[V](fi) traces in Fig. 0in terms of the peculiar but 
regular vq-vs-H pattern sketched in Fig. [|J] 



III. HALL BAR GEOMETRY 

A. Boundary conditions and kernels 

We now consider a 2DEG with lateral confinement in 
the x direction and translation invariance in the y di- 
rection, i.e., an idealized Hall bar geometry. To study 
boundary effects on the screening properties, we will ap- 
ply an additional periodic external modulation poten- 
tial in x direction. We will consider two different sets 
of boundary conditions, which lead to slightly different 
confinement potentials. 



In-plane gates 



Following Refs. ITollnllTlT^ we first assume that all 
charges reside in the plane z = 0, and that the halfplancs 
z = 0, x < —d and z — 0, x > d are kept at constant 
electrostatic potential, V(x, y, z = 0) = for \x\ > d (in- 
plane gates) i 11 ! 13 Then the electrostatics can be solved 
using the theory of complex functions, and the kernel in 
Eq. O, with —xi = x r — d, is obtained as^ 3 - 



K l{ (x,t) 



In 



v/(d 2 -x 2 )(d 2 -t 2 ) + d 2 



fx 



(x - t)d 



(28) 



Positive background charges of the 2D charge density eno 
between the in-plane gates will produce the confinement 
potential (written as potential energy of an electron) 



Vb g (.x) = -Eoy/1 - (x/d) 2 , E Q = 2Tre 2 n d/R, (29) 

which can be calculated from Eq. using the kernel 
and replacing n e i(x') by — no- 
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FIG. 8: Confinement potential Vb g (x,y,z — 0)/Eo due to a 
homogeneous plane of charge density eno at distance z. The 
dash-dotted line is obtained from model 128H with z = 0, the 
other lines are for model (1301 . 



2. Perpendicular gates 

Another simple set of boundary conditions is obtained 
assuming the 2DEG to be laterally confined by two 
equipotential planes located at x = ±d parallel to the 
y-z plane, V(x — ±d,y,z) = 0. This is a reasonable 
model for a free standing mesa-etched Hall bar with free 
or metallized surfaces at x — id, which accommodate 
a large number of (partially occupied) surface states. 
The electrostatics with these boundary conditions is well 
known^ In our notation it is expressed by the kernel 



K x{ x,t) = -J C ° S l^ X + t)+ ^ 



(30) 



for 7^0. Inserting this with 7 = sinh(7rz/4<i) into 
Eq. ||TJ, where — xi — x r — d, yields the electrostatic 
potential Vh{x, y, z) due to the 2DEG at a position sep- 
arated by the distance \z\ from the plane of the 2DEG. 
Correspondingly, we can use this to calculate the con- 
finement potential produced in the plane of the 2DEG 
by a plane, positive background charge at a distance 
z from the 2DEG. Typical confinement potentials are 
shown in Fig. |HJ For 7 = the potential minimum is 
V hs (0,y,0)/E = -8G/71- 2 = -0.74246, with Catalan's 
constant G = 0.915965594. 

The positive background charge density and the sam- 
ple width define the characteristic energy Eq, Eq. H29fl . 
Measuring energies in units of Eq, lengths in units of 
d, and density of states in units of Dq, we obtain from 
Eq. the dimensionless electron density n(x/d) — 
n e \(x)/ EqDq, so that Eq. Q assumes a dimensionless 
form with the prefactor l/a con f, where 



a con f = na* B /2d 



(31) 



measures the relative strength of the Coulomb interac- 
tion, similar to a of Eq. 10. We will usually assume 
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FIG. 9: Some consistent density profiles (upper panel) 
and potentials (lower panel) calculated for the in-plane-gates 
model (1281 and the perpendicular-gates model (1301 . respec- 
tively. The depletion length is chosen as d/5 for the short- 
dashed curves and as d/10 else. Thin horizontal lines indicate 
the corresponding electrochemical potentials, a — na B /2d, 
T = 0, B = 0. 

aconf = 0.01, i.e., for GaAs, a sample width 2d ~ 3^m, 
since this allows us to calculate density profiles with 
clearly visible incompressible strips on a mesh of rela- 
tively few (~ 500) points across the sample. For much 
larger d, we would need a much finer mesh, i.e., more am- 
bitious numerics, and the incompressible strips would be 
hardly visible on that scale, although the physics would 
not change qualitatively. 

Figure [3] shows some density and potential profiles 
obtained for the two sets of boundary conditions in 
the limit of zero temperature and magnetic field, where 
n cl (x)/n = (Tr/a con{ )n(x)0(^(x))/E , with fi(x) = fi* - 
V(x). Apparently the density profiles are very similar, if 
we assume the same depletion length, the same sample 
width, and vanishing spacer between 2DEG and back- 
ground charges (i.e., z = 0) in both cases. 



B. Unmodulated system in a magnetic field 

In the ideal homogeneous 2DEG at low tempera- 
tures, the chemical potential as a function of the mag- 
netic field exhibits the well known saw-tooth behavior, 
Eq. (|27|l . With decreasing B it follows a Landau energy 
(n+ l/2)frui c until the filling factor v = 2Ef jhw c reaches 
the value 2(n + 1), and then it jumps to the next higher 
LL. In the confined system, the self-consistently calcu- 
lated "chemical potential" ^t(0) = fi*— V(0) in the center, 
x = 0, shows the same behavior, as is seen in Fig. IKJL . 
where (j,(x = 0;i3,T) in units of fig = fj,(x — 0;0,0) is 
plotted as function of £1 = fkt) c / jiQ. 
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FIG. 10: a) Filling factor v{x) and b) potential V(x) for 
the values of Q = HiOc/^o given in the legend and indicated 
by open circles in c), which shows the "chemical potential" 
fj, — (i* — V(0). Calculated with model 1281 for parameter 
values ira* B /2d = 0.01, no/E = 0.00284, k B T/E = 2 x 10 -5 . 

However, in contrast to the chemical potential oscil- 
lations in a homogeneous 2DEG, the corresponding os- 
cillations in the confined system are realized by strong 
spatial variations of the electrostatic potential in the in- 
terior of the sample. This is demonstrated in Fig. 110b . 
which shows the self-consistent total potential in the in- 
terior of the sample for the four values of fi indicated 
by open circles in Fig. HOh. Figure Unk shows the corre- 
sponding density profiles, normalized as local filling fac- 
tor, v(x) = 2Trl 2 n cl (x). For v(x) < v(0) « 2/fi < 2 
the LL n = is pinned to the electrochemical potential 
fi* nearly everywhere in the 2DEG. For v(0) > 2 the 
LL n = 1 must be partially populated in the center of 
the sample. This forces the potential to develop a lo- 
cal minimum near the center, with a decrease of V(0) 
by an amount ~ fiw Cl so that a compressible strip starts 
to develop in the center. With further increasing ^(0), 
this central compressible strip becomes broader and the 
adjacent incompressible strips, together with the related 
potential steps, move towards the sample edges. Similar 
drastic changes of the potential distribution are found 
near all jumps of the chemical potential. Thus, pinning 
and screening lead already to drastic effects in the con- 
fined 2DEG even in the absence of any additional poten- 
tial modulation. 



C. Confined system with modulation 

We now add a symmetric external modulation poten- 
tial to the confinement potential and investigate how this 
affects the self-consistent potential. We take V^xtO^) = 
Vq cos(2.57ra;/d) which is in accord with our general 
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boundary conditions and exhibits just one full oscillation 
period in the interior of the sample, so that we can expect 
similar screening effects as in a homogeneous unbounded 
system, and possibly some effects of the nearby sample 
edges. The period of this modulation is a = 2d/2.5, so 
that the choice a CO nf = 0.01 implies a = 1/40 [see Eq. (JjJJ] 
and the results can be compared immediately with the 
previous one for the unbounded 2DEG. For this compar- 
ison it will be important whether the potential of the 
unmodulated system has a strong variation in the center 
region or not, i.e., whether the filling factor ^(0) in the 
center is slightly larger than an even integer or not. 




FIG. 11: a) Filling factor v(x) and b) screened poten- 
tial AV{x;Vo) for the confinement models Ij28^ (solid lines) 
and 1301 1 (dash-dotted lines), and for the unconfined 2DEG 
(dashed lines). Thin lines show v(x) without modulation, 
with v(0) = 2.8574 for all models. The screened potential for 
the unbounded 2DEG is shifted by a constant, and actually 
oscillates symmetrically around zero. (a CO nf = 0.01, Ho/Eo = 
0.002842, hujc/no = 0.7, V /no = 24.63, k B T/[i = 0.007) 



1. Weak boundary effects on screening 

To describe the screening of an external modulation 
potential, it seems natural to calculate the difference 
AV(x; Vq) = V(x; Vq) — V(x; 0) of the self-consistent po- 
tentials with and without the modulation. If V(x; 0) is 
flat in the interior of the sample, we expect that screen- 
ing is very similar to that in an unconfined system and 
that AV(x; Vq) contains essentially the same information 
as V(x;Vo), apart from an unimportant constant offset. 
This is indeed true if z/(0) is not closely above an even 
integer. As an example, we compare in Fig. m numerical 
results for the two confinement models and for the uncon- 
fined 2DEG. The results for the density of the confined 
2DEG differ only slightly in the edge regions. In the inte- 
rior, the filling factors v{x) are practically the same and 



agree well with that of the unconfined 2DEG with the 
same modulation potential. Also the screened potentials 
are equivalent and differ only by a constant offset, which 
results from the asymmetry of the density modulation 
with respect to the unmodulated electron density profile. 



2. Strong confinement effect on screening 

Things become more complicated, if already without 
additional external modulation the potential near the 
center of the Hall bar varies strongly, as happens for 
v(Q) = 2n + v n with < v n <C 1. Then, for small 
modulation amplitude (Vq <C hui c ), the self-consistent 
potential V(x;Vq) follows V(x;0), with a minimum at 
x = 0, and only the difference AV(x; Vq) reminds of an 
oscillatory potential with the phase of the external mod- 
ulation, see solid lines in Fig. For stronger modu- 
lation (Vq < Sw c ), V(x;Vq) develops a local maximum 
at x — and the total variation of V(x; Vq) in the cen- 
ter region |x| < d/2 is of the order of Vq < hoj c see 
Fig. 112b. The variation of AV(x; Vo) is now, however, by 
an amount of Hlo c larger. In this small- Vq regime screen- 
ing is rather poor and very nonlinear. As Vq increases 
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FIG. 12: a) Selfconsistent potential V(x; Vo) and b) screened 
potential AV(x; Vo), for several Vo. For clarity, b) shows 10 • 
AV(x; Vo) for the weakest modulation Vo = 6 • 10 -4 -Eo- For 
Vo = [thin solid line in a)] u(0) = 2.03. (a conf = 0.01, 
Mo/So = 0.002842, hw c /fi = 0.985, k B T/fi = 0.01) 

further, the variance Var[F](F ) = V(0; Vo) - V(a/2; V ) 
(note that a/2 — 0.4 d) approaches the plateau value fkt! c , 
and then behaves as a function of Vo just as for the un- 
confined 2DEG. The variance of the "screened potential" 
AV(x; Vq), on the other hand, will be by about fkj c larger 
in the plateau region. 

To summarize, if we neglect the relatively narrow Vq 
interval between the plateaus, we find that the variance of 
the self-consistent potential V(x; Vq) as a function of Vq 
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shows the same behavior as for the unconfined 2DEG. For 
most magnetic field values, the variance of the "screened 
potential" AV(x; Vq) also shows the same characteristics. 
Only if the filling factor i/(0) in the center of the unmod- 
ulated confined 2DEG is slightly larger than an even in- 
teger, the spatial variation of the self-consistent potential 
V(x; 0) of the unmodulated system causes the variance 
of AV(x; Vq) to be about huj c larger than the variance of 
V(x; Vb). If we plot the variance of AV(x; Vo) at fixed Vo 
as a function of f2 = Hu) c /fio, we get sawtooth-like traces 
as in Fig. [7| however with additional spikes of height 
at fi < 1/k for integer k. 

IV. SUMMARY 

We have investigated the screening of a harmonic ex- 
ternal potential by an unconfined two-dimensional elec- 
tron gas as well as by confined 2DEGs in a simplified 
Hall geometry, in strong perpendicular magnetic fields 
and at low temperatures. Our numerical results within 
the self-consistent Thomas-Fermi- Poisson approach show 
that screening is very nonlinear and dominated by the 
phenomenon of pinning of Landau levels to the electro- 
chemical potential, which leads to compressible regions 
with position-dependent electron density, where this pin- 
ning takes place, and to incompressible regions of con- 
stant density and position-dependent electrostatic po- 
tential in between. At fixed magnetic field, the total 
variation ("variance" Var[V]) of the self-consistently cal- 
culated potential energy increases with increasing mod- 
ulation amplitude Vo in a step-like fashion, exhibiting 
plateaus, where the value of Var[V] is close to an inte- 
ger multiple of the cyclotron energy and shows a weak 
linear increase with Vo, with a slope proportional to the 
temperature. The corresponding modulation of the elec- 
tron density is, in contrast to the potential, not strongly 
affected by the magnetic B. The occurrence of incom- 



pressible strips leads to local modifications, but the over- 
all density profile is roughly the same as for B = 0, as 
has already been emphasized by Chklovskii et aliH. Ex- 
ploiting this observation together with the pinning phe- 
nomenon and the relations between density modulation 
and external and screened potential valid in the linear 
screening regime, we were able to derive simple analyt- 
ical expressions for step heights and plateau widths of 
the Var[V~]-vs-Vb curves for arbitrary B and T = 0. This 
simple analytical description of nonlinear screening in an 
unconfined 2DEG is summarized in Fig. [U and allows 
also an easy understanding of the complicated traces ob- 
tained while plotting Var[V~] as a function of B at fixed 
Vo (see Fig. E). 

Finally we have investigated the corresponding screen- 
ing properties of a confined 2DEG in a simplified Hall 
geometry, for two different types of boundary conditions, 
which lead to different confinement potentials, but nearly 
identical density profiles, apart from slight deviations in 
the edge regions. Considering the effect of an external 
modulation potential V ext (x) = Vo cos(2.57ra;/a) in the 
interior of the sample (more than about a/2 from the 
edges), we find essentially the same properties as for the 
unbounded 2DEG. Care must however be taken if in the 
center of the unmodulated system a new Landau level 
starts to be occupied, since then the self-consistent po- 
tential varies strongly in the center region. This is an 
interesting confinement effect, but it can be easily elimi- 
nated from the discussion of screening if the modulation 
amplitude Vo is large enough. 
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